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Abstract 

In this short note we show how virtual arbitrage opportunities can be mod- 
elled and included in the standard derivative pricing without changing the general 
framework. 

Whatever people say about the drawbacks of the Black-Scholes (BS) approach |l|] 
to derivative pricing, it is a standard method and almost any pricing and hedging 
software in financial institutions is based on it. Practitioners have got used to BS- 
like partial differential equations, martingales and other related mathematical animals. 
Both analytical and numerical methods are well developed and it is hardly surprising 
that practitioners are rather reluctant to "buy" complicated new theories. That is why 
it is interesting to see how some limitations of BS analysis can be overcome in the same 
mathematical framework without disturbing the foundations. 

One way to improve BS is to use a more realistic price process instead of the 
geometrical random walk. The most popular alternatives are ARCH-GARCH models 
where the volatility of the return is assumed to be stochastic. Although they are 
a better approximation for the price process the description is far from perfect 0. 
Another line of attack is the no-arbitrage constraint. Some empirical studies have 
demonstrated existence of the short lived arbitrage opportunities @, 0. In this paper we 
show how to generalize BS analysis to the case where the virtual arbitrage opportunities 
exist. 

In Ref Ilinski and Stepanenko proposed BS-based model to account for virtual 
arbitrage. The model substitutes the constant interest rate Tq as a rate of return on 
a riskless portfolio by some stochastic process tq + x{t). In the limit of fast enough 

*E-mail: kni@th.ph.bham.ac.uk 



market reaction to the arbitrage opportunities (see [§ for details) the model is reduced 
to the Ornstein-Uhlenbeck process for the variable x{t): 

^ + Aa; = vit) (1) 

where rj(t) is the white noise stochastic process 

{Vit))=0, {r,{t)r,{t'))=T?5{t-t'). 

The parameter A characterises how fast the market reacts to an arbitrage opportunity 
while the parameter E^ defines how often such opportunities appear and how profitable 
they are. In the limit S^ ^ or A — > oo we recover the results of the BS analysis. 

As was already mentioned, real price processes are not geometrical random walks. 
However, if we assume they are, the parameters S^ and A can be estimated from the 
market data as 
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with the Black-Scholes operator £bs defined as usual: 

Here 5* is an underlying asset price and a is its volatility. The (historic) market values of 
the parameters allow us to estimate how appropriate the no-arbitrage approximation is, 
i.e. when it is possible to neglect the correction terms come from the virtual arbitrage. 
To first order with respect to the parameter ^ the correction to the BS equation was 
found in Ref [§] together with explicit formulas for the call and put vanilla options. 
However it is clear that in the case ^ > 1 this expansion does not make much sense 
and we have to find a way to sum over all relevant terms. 

How real is this possibility? Let us assume that the arbitrage on the riskless portfolio 
V —S^ appears every day and is on average about one per cent (in absolute value). We 
assume also that the arbitrage is washed out during the day and there is no interference 
between different arbitrage opportunities. Then the parameter A is about 260, the 
arbitrage produces the return of about 260 per cent per year and the parameter ^ 
is about 6.76 = (2.6)^ >> 1. In this case the no-arbitrage BS pricing formulas are 
too crude an approximation. If however the arbitrage appears only once a week then 
^ ~ 0.25 and the parameter can be used as an expansion parameter. These examples 
show that accounting for the virtual arbitrage can sometimes be very appropriate. It 
is worth mentioning that here we have neglected the bid-ask spread and transaction 
costs for illustrative purposes. In reality both effects have to be included. 



It is not difficult to derive an analogue of the BS equation for the derivative price 
U{S,t) averaged over virtual arbitrage process (|I]): 



U{S,t)= dxu{x,S,t) (2) 

where the function u{x, S, t) is the solution of the problem: 
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u{x,S,t)\t=T = Six) ■ Payoff (S) . 

A concise informal derivation of the equation is in the Appendix. It is straightforward 
to check that U{S, t) converges to the solution of the BS equation as S^ or 1/A tend to 
zero. It has to be used for both pricing and hedging, allows standard numerical methods 
to be used in its calculation and can serve as a basis for further generalizations. 
Two such generalizations are immediately obvious: 

1. It is possible to use another form of the stabilizing market reaction, if it fits the 
observable data better. In this case the process will be defined as 
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with the same white noise r]{t) 

{V{t))=0, {r,(tW))=Y?5{t-t'), 

which leads to the following problem for the function u{x, S, t): 

u{x, S, t)\t=T = 5{x) ■ Payoff (S) . 
The average price lf{S,t) then is found as before from Eqn(0). 

2. In a similar manner it is possible to improve the stochastic volatility models |^, 
1^, H, the jump diffusion model and others. For example, in the case of the 
price process S with stochastic volatility a: 

dS = (j)Sdt + aSdw , V = a\ dV = fiVdt + ^Vdz 

such that the Wiener processes have correlation p and the volatility process has 
zero systematic risk 0, the pricing equation takes the form: 
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with the corresponding final condition: 

m(x, S, t)\t=T = S{x) ■ Payoff (S) . 
The average value of the price is given by Eqn(|^). 



It is possible to combine these generalizations or expand the treatment to other 
BS-like equations. In any case the choice of one variant or another should involve 
extensive empirical study. It is important to emphasize however that all the pricing 
ideology, technology of calculations and numerical methods remain essentially the same 
so that the corresponding software can be easily adapted for the presence of the virtual 
arbitrage. 

Lastly we note an apparent analogy between the virtual arbitrage return and a 
stochastic interest rate with mean reversion (I am grateful to Neil Johnson and Stefan 
Thurner for this comment). However, as it can be seen from the numeric values, the 
interest rate time-scale is much larger while the characteristic fluctuation size is con- 
siderably smaller than the virtual arbitrage counterparts. This makes the corrections 
valuable even for short living derivatives when the stochasticity of the interest rate 
can be neglected. Another difference is more theoretically important. In the models 
with stochastic interest rate the risk can be actually hedged using bonds with various 
maturities. In the case of the virtual arbitrage this risk is intrinsic and cannot be 
hedged. In other words there is no tradable instrument involving the same risk which 
can be used for pricing and hedging. Even the whole concept of pricing and hedging 
using duplication becomes somewhat obscure. 
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Appendix 

In this appendix we give an informal derivation of the main equation (^ using the 
formahsm of functional integrals. 

We start with the Black-Scholes equation with the rate of return on the riskless 
portfolio V — S^ given by ro + x: 

^ + ^-^ + (ro + x)5— -(ro + x)I/ = 0, 

V{S,t)\t=T = S{S-S') . 

To simplify the calculation we change variables as r = T — t and y = In S which casts 
the previous equations in the form: 
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The solution of the problem can be expressed as the following functional integral |T0 



V(y,T) = - f DqD^ j;{ipq-4p^-{ro+x(r))+i{ro+.{r)-^)p)dr 
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with the boundary conditions g(0) = y', q{T) = y. The functional integral form of 
the solution is extremely convenient for the purpose of averaging over trajectories x{t) 
since it presents the dependence in explicit form. We, however, will use another trick. 
Let us first of all average the previous expression over realisations of the stochastic 
process x(r) with fixed ends, i.e. when x(r) = x and x(0) = (since there is no 
arbitrage at the expiration date of the contract and later). The probabilistic weight 
of a particular trajectory x(r) for the Ornstein-Uhlenbeck process (|I|) is given by the 
expression 
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which leads to the following result for the conditional average value of the contract 
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?(0) = y , g(r) = y' , x(0) = , x{t) = x . 

Now, instead of evaluating these integrals, we find a partial differential equation for 
Vix,y,T) 11|: 

dV a^d^V , a'^.dV , ,- Y? d'^V 3xV 



with the initial conditions 
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Returning back to the initial variables t and S we obtain the problem (^) 

Integration of the solution over x (to get the unconditional average) and the convolution 
with the final payoff complete the consideration. 

In a similar way various generalizations of the equations can be derived. It is clear 
that the functional integral is a very convenient tool for such kind of manipulations. 
We would like to note however that the functional ingral formalism is full of subtitles 
which we have not emphasized here. It serves for fast derivation but not proving the 
results. The proof can be obtained in the routine framework of stochastic calculus. 



